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Abstract 

We present the calculations of the complete next-to-leading order(NLO) QCD corrections to the 
inclusive total cross sections for the Kaluza-Klein(KK) graviton and photon associated production 
process pp — > jGkk + X in the large extra dimensions model at the LHC. We show that the NLO 
QCD corrections in general enhance the total cross sections and reduce the dependence of the total 
cross sections on the factorization and renormalization scales. When jet veto is considered, the 
NLO corrections reduce the total cross sections. We also calculate some important differential cross 
sections for this process at NLO: the missing transverse momentum distribution, the transverse 
momentum distribution and the pseudorapidity distribution of photon. 
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I. INTRODUCTION 



The idea that the extra dimensions theory can appear at the 



eV scale, well below 



lMlO|| and promises rich 



the Planck scale Mp ~ 1.2 x 10 19 GeV, was proposed in the 1990s 
phenomenology at the TeV scale. Now the search for extra dimensions is one of the important 
tasks for the LHC. 

Among various extra dimensions models the large extra dimensions(LED) model in- 
troduced in Ref.[l-3] is the first TeV scale gravity theory and has been extensively 



studied 
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In this model space-time has 4 + 5 dimensions and the standard model(SM) 
particles reside in the usual 3 + 1-dimensional SM brane and can not propagate in the extra 
5-dimensional space, which is assumed to be compacted on a torus with a common radius 
R, while gravity can propagate in the whole 4 + 5 dimensional world. From the view of our 
4-dimensional world, there exists infinitely many Kaluza-Klein(KK) modes of gravitons with 
mass \k\/R interacting with SM particles, where k 2 = Ylt=i^i w hh ki being the integer. 
The 4-dimensional Planck scale is no longer the fundamental scale, but an effective scale in 
the 4-dimensional world, and is related to a fundamental scale Mp ~TeV by the Newtonian 
law of gravitation in 4 + 5 dimensions 1, [l^] 

M P = Mp/V^ = R S/2 M 1 D +S/2 , (1) 

where Mp is the reduced Planck mass. According to Eq.flT]), deviations from the usual 
Newtonian gravitational law appear at R ~ 10^~ 19 meters, which is not a conflict with 
the current gravitational experiments 13J] once 8 > 2. Before further results of terrestrial 
gravitational experiments appear, one may resort to colliders to find signals of this model. 
Although, in this model the couplings of gravitons to SM particles are suppressed by 1/Mp 
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121 ]. the summation of the production of large numbers of KK modes with arbitrary 



mass smaller than Mp, may compensate for this suppression and lead to observable effects. 
There are two ways to probe such effects at collider s: g raviton emission and virtual graviton 



exchange, which have been investigated in Ref. ll.fl^fl^. As shown in Ref. 12], because of 
the suppression of the couplings of gravitons to SM particles, the decays of gravitons to SM 
particles have a small probability to occur before they propagate into the extra £- dimensional 
space, which means the gravitons behave like massive, stable, and noninteracting particles 
once they are produced. Thus, the signal for graviton and photon associated production at 
the LHC is a single photon plus missing energy. Since the electromagnetic coupling is small 



and the qq luminosity is lower than for gg at the LHC the rate for this process is much 
smaller than for jet and graviton associated production. But the photon signal would be a 
clean signature; and in case of discovery in jet plus graviton events, the photon plus graviton 
signal would provide a useful independent test [12]. Only leading-order (LO) calculations and 
analysis of the process were performed in Ref.[12j. Since LO cross sections for processes at 
hadron colliders suffer from large uncertainties due to the choices of the renormalization scale 
(fi r ) and factorization scale (///) and higher order QCD corrections are generally large and 



can improve the scale uncertainties at hadron colliders. Several works |15H23| have performed 
next-to-leading order(NLO) QCD corrections in extra dimensions models. In this paper we 
present the complete calculations of NLO QCD corrections to this process which improve 
the theoretical predictions. 

This paper is organized as follows. In SecJTT] we show the analytic results of the LO 
calculations and define the notation. In Sec JIIII we present the details of the calculations of 
both the virtual and real parts of the NLO QCD corrections. In Sec JIVI we give the numerical 
predictions for inclusive and differential cross sections at the LHC. We close this paper with 
a brief conclusion. For completeness, the relevant Feynman rules are collected in Appendix 
A and the lengthy analytic expressions of the results of our calculations are summarized in 
Appendix B. 

II. LEADING-ORDER CALCULATIONS 

The KK gravitons with different masses can be produced at colliders only if kinematically 
allowed. Contributions from the different KK modes then must be summed up. Since the 
KK graviton mass separation of 0(1/ R) is much smaller than all the other physical scales 
involved, we can replace the discrete summation of different KK modes by a continuous 
integration. In general, the differential cross section for graviton production can be expressed 



as 



12| 



d 2 a M 2 P * x da m 



dt dm M^ +5 dt 

with 

2tt 5 / 2 

Ss - 1 = TW2? (3) 

where Ss-i is the surface of a unit-radius sphere in 5 dimensions and da m /dt is the differential 
cross section for producing a single KK graviton of mass m. Throughout this work we 
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FIG. 1: Leading-order Feynman diagrams for qq — > 'jGkk- 

perform the integration on m from O.IMjj to Mjj. 

The leading-order Feynman diagrams for the graviton and photon associated production 
process q(p-[)q(p2) — > ''iip^GKKipi) are shown in Fig. CD The related Feynman rules are 



given in Ref. 
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in Ref. 



121 ] and are collected in Appendix A. The LO amplitudes have been given 



111 ] so we need only to show the amplitudes squared here. In the LED model the 



spin sum over the polarization tensors of the graviton is 

5 



E s s* _ p 



with 



n 



(4) 



(5) 



where the dots represent terms proportional to the graviton momentum and p^T^ = 0, 
giving no contribution to the amplitude. We performed our calculations in n = 4 — 2e 
dimensions. 



The LO partonic cross section is then 



2s 



B\2 



(6) 



with 



2 ,.2 



B\2 



eQ k 



[{s + At)m 6 - Qt(s + 2t)m A + (s 3 + Qts 2 + 18t 2 s + lQt 3 ) 



2Astu 

-At (s 3 + 3ts 2 + 4t 2 s + 2t 3 )] 



m 



where k = y2/Mp, s, t, and u are the Mandelstam variables defined as 

8 = (Pi + P2) 2 , t=(px- p 3 ) 2 , u = (pi - Pa) 2 , 



(7) 



(8) 



and means that the colors and spins of the outgoing particles have been summed over 
and the co.ors and spins of the incoming paries have been averaged over. Equation© is 
a coincidence with the result shown in Ref . [12[ • 

The LO total cross section can be obtained by convoluting the partonic cross sections 
with the parton distribution functions(PDF) G Qyq in the protons: 

°"m = J dx 1 dx 2 [G q/p (x 1 , fi f )G 9 / p (x 2 , n/) +G q / p (x 2 ,Hf)G 9 / p (x 1 ,fif)]<T%. (9) 
Here ///is the factorization scale. 



III. NEXT-TO-LEADING ORDER CALCULATIONS 

The NLO corrections to the associated production of a graviton and a photon can be sepa- 
rated into the virtual corrections arising from loop diagrams of colored particles and the real 
corrections arising from the radiation of a real gluon or a massless (anti) quark. We carried 
out the calculations in the 't Hooft-Feynman gauge and used dimensional regularizationj^j] 
in n = 4 — 2e dimensions to regulate all the ultraviolet (UV), soft and collinear divergences. 
We performed two independent calculations for both the analytical and numerical results 
for cross checking, and the results of the two groups agree with each other. 



A. Virtual corrections 



The Feynman diagrams for the virtual corrections to qq — > ^Gkk are shown in FigfS] and 
Figj3] They consist of self-energy, vertex, triangle and box diagrams. In order to remove 
the UV divergence we adopt the on-shell renormalization scheme 25N28]. 

We denote the bare and renormalized quark wave functions by tp q o and ip q , respectively. 
The renormalization constant 6Z q is then defined by 



i> q0 = (l + 6Z q ) 1 / 2 i> q 



(10) 



Calculating the quark self-energy diagram we obtain the explicit expression for 5Z q : 

SZ q = ^C F (- - — ). 
4:71 e e uv 

Here Cf = f while 1/e and I /euv represent infrared(IR) and UV divergences, respectively. 



(11) 



After renormalization the UV divergences in the virtual corrections are removed leaving 
the IR divergences and the finite terms. The 0(a s ) virtual corrections to the partonic total 
cross section can then be expressed as 

with 

A% = -C F , A\ = -\c F . (13) 

Here S represents finite terms in the virtual corrections and the explicit expressions are 
given in Appendix B. The cancellation of IR divergent terms 1/e 2 and 1/e will be discussed 
in detail below. 



o:. s 



r(i 



27rr(l-2e) 



f2A\ 



2 _|_ ^^l 



|M 



B\2 



+ —S 

2vr 



(12) 



B. Real gluon emission 

The Feynman diagrams for the real gluon emission process q{pi)q{p2) — > i{Pz)Gkk{P's) + 
g{p$) are shown in Fig. HI 

The phase space integration for the real gluon emission will produce soft and collinear 
singularities which can be conventionally isolated by slicing the phase space into differ- 
ent regions using suitable cutoffs. In this paper we use the two-cutoff phase space slicing 



method 29|], which introduces two arbitrary small cutoffs; i.e., a soft cutoff 5 S and a collinear 
cutoff 8 C , to divide the three-body phase space into three regions. 

First, the phase space is separated into two regions by the soft cutoff 8 S , according to 
whether the energy of the emitted gluon is soft; i.e., E 5 < 5 sy /s/2, or hard; i.e., E 5 > 5 sy /s/2. 
Then the parton level real cross section can be written as 

= + ( 14 ) 

where and a% are the contributions from the soft and hard regions, respectively, 
contains all the soft divergences, which can explicitly be obtained after the integration over 
the phase space of the emitted gluon. Next, in order to isolate the remaining collinear 
divergences from a^, the collinear cutoff 5 C is introduced to further split the hard gluon 
phase space into two regions, according to whether the Mandelstam variables = (pi— P5) 2 , 
with i = 1,2, satisfy the collinear condition — 5 c s < t i5 < or not. We then have 

(J m= a ra + °m , ( 15 ) 




FIG. 2: One-loop virtual diagrams including vertex and self-energy corrections to qq — > 'jGkk- 
Each brown vertex is UV divergence free. 

where the hard collinear part contains the collinear divergences, which also can explicitly 
be obtained after the integration over the phase space of the emitted gluon. The hard 
noncollinear part is finite and can be numerically computed using standard Monte 



Carlo integration techniques 30] and can be written in the form 

d^=^\M^dT,. (16) 

Here dT 3 is the hard noncollinear region of the three-body phase space. 

In the next two subsections we will discuss in greater detail the soft and hard collinear 
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FIG. 3: Box and triangle diagrams for qq — > •jGkk- The UV divergences cancel among the five 
diagrams. 

gluon emission. 

1. Soft gluon emission 

In the limit that the energy of the emitted gluon becomes small, i.e. E 5 < 5 sA /s/2, 
the amplitude squared J2\M(qq — > ■jGkk + fiO| 2 can be factorized into the Born amplitude 
squared times an eikonal factor <3> eik : 

Y^\M{qq -> iGkr + g )\ 2 ^ (47r« s/ i^)^|M B | 2 $ eik , (17) 
where the eikonal factor is given by 

$cik = C F -. £ -. (18) 

(Pi •PB)[P3-P5) 

Moreover, the three-body phase space in the soft limit can also be factorized: 

dT 3 (qq -> -yG KK + g)^% dT 2 (qq -> ^G KK )dS. (19) 

n 

Here is the integration over the phase space of the soft gluon and is given by|29j 

1 r-Ssy/s/2 

dS = tttttt^t / dE 5 Et^dQ 2 ^ 2e . (20) 



2(2vr; 



o 
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FIG. 4: Feynman diagrams for qq — > 'jGkk + 9- 



The parton level cross section in the soft region can then be expressed as 



a s m = (47ra s/ i r 2e ) J dT 2 J2\M B \ 2 J dS$ cik . 



(21) 



Using the approach in Ref.|29|, after integration over the soft gluon phase space, Eq.pT 
becomes 



° m ° m 



a s T(l — e) / 47r/i 



2ti T(l - 2e) 



,2\ e 



-4 + — + A 



with 



A s 2 = 2C F , A{ = -AC F log S s , A s = AC F log 2 S s 



(22) 



(23) 



2. Hard collinear gluon emission 



In the hard collinear region, E 5 > 5 s y/s/2 and — 5 c s < t& < 0, the emitted hard gluon is 



collinear to one of the incoming partons. As a consequence of the factorization theorem [31 



32] the matrix element squared for qq — > , jGkk + <? can be factorized into the product of the 
Born amplitude squared and the Altarelli-Parisi splitting function for q(q) — > q [33N37I| : 
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that is, 



Y)M{qq -+ !G KK + g)f C °^T (4vra s ^)^|M B | 2 (-^^ + ^f^ ) > ( 24 ) 

where z denotes the fraction of the momentum of the incoming parton carried by q(q) with 
the emitted gluon taking a fraction (1 — z). Pij(z, e) are the unregulated splitting functions 
in n = 4 — 2e dimensions for < z < 1 which can be related to the usual Altarelli-Parisi 
splitting kernels [33] as follows: Pij(z,e) = Pij(z) + eP[Az\ Explicitly 



P m {z) = P m {z) = C F ^^ + C F ~6(1 - z), (25) 
P' qq {z) = P'-^z) = -CV(1 -z) + C F U{1 - z). (26) 

Moreover, the three-body phase space can also be factorized in the collinear limit and, for 
example, in the limit — 5 c s < ti§ < it has the following form [29]: 

dT 3 (qq -> jGkk + fiO c °^T r dT 2 (q(q) -> jGkk] = zs) -dzdt 15 [-(l - z)t i5 ]~ e . 

l07T z l (1 — 6) 

(27) 

Here the two-body phase space should be evaluated at a squared parton-parton energy of 



zs. Thus the three-body cross section in the hard collinear region is given by 29j 

'a s r(l-e) /4tt/i 2X 6 



dal c = da, 



~~) S c 6 [Pqq(z,e)G q /p(x 1 /z)Gg/ p (x 2 ) 



2nT(l -2e) 

dz (\ - z\ ~ e 

Pqq(z, e)Gq/ p (x 1 /z)G q /p(x 2 ) + (x\ <-> x 2 )] — ( j dx x dx 2 (2£ 



where G q ( q y p (x) is the bare PDF. 



C. Massless (anti) quark emission 

In addition to real gluon emission a second set of real emission corrections to the inclusive 
cross section for pp — > •jGkk at NLO involves the processes with an additional massless q(q) 
in the final state: 

q{q)g^lG KK + q(q). (29) 

The relevant Feynman diagrams are shown in Fig. |5j The diagrams for q emission are 
similar and are omitted here. 
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FIG. 5: Feynman diagrams for qg — > jGkk + Q- 



Since the contributions from real massless q(q) emission contain initial state collinear 
singularities we need to use the two cutoff phase space slicing method [29] to isolate these 
collinear divergences. But we only split the phase space into two regions because there are 



no soft divergences. Consequently, using the approach in Ref. 



291 ]. the cross sections for the 



processes with an additional massless q(q) in the final state can be expressed as 



da. 



add 



^2 ^Zi a P iGkk + q{q))[G a / p {x 1 )Gp/ p {x 2 ) + (x 1 x 2 )]dx 1 dx 2 

(a=g,/3=q,q) 



+da. 



B 



a s T(l - e) 



47T/1 



2\ e ' 



2?r TCI — 2e) \ s 

+Pqg{z,e)G q / p (x 1 )G g / p (x 2 /z) + (x 1 x 2 )] 



~) s c 6 [P qg (z,e)G g/p (x 1 /z)G q/p (x 2 ) 
dz ( 1 



dx\dx 2 



(30) 



where 



l\z 2 + (l 



-z(l-z). 



(31) 



The a^ term in Eq. fl30|) represents the noncollinear cross sections for the q(q)g initiated 
processes which can be written in the form 



noncollinear 



lG KK + q(q))\dT 



(32) 



11 



where dT^ is the three-body phase space in the noncollinear region. The other terms in Eq. 
fl30|) are the collinear singular cross sections. 



D. Mass factorization 



After adding the renormalized virtual corrections and the real corrections, the parton level 



cross sections still contain collinear divergences which can 
the PDFs at NLO, generally called mass factorization 38 



3e absorbed into a redefinition of 
. This procedure, in practice, 



means that first we convolute the partonic cross section with the bare PDF G a / P {x) and 
then use the renormalized PDF G a / P (x, /i/) to replace G Q / P (x). In the MS convention the 
scale- dependent PDF G a / P (x, ///) is given by [291 ] 



G a/p {x,fi f ) = G a/p {x) + (--) 



a s T(l - e) 
27rlYl - 2e) 



x 



47T/X 



- 1 dz 

x / —P a p{z)Gp/ p {x/z). 



(33) 



This replacement will produce a collinear singular counterterm which is then combined with 
the hard collinear contributions to give Ref. [29J] the 0(a s ) expression for the remaining 
collinear contribution: 

~a s r(l-e) /4vr/i^ e T 



da™ 11 = da% 



2tt r(l - 2e) v « 
Af(a -> ag) 



{G q / p (Xi, /J,f)Gq/p(x 2 , Hf) + G q / p (Xi, flf)Gq/ p (x 2 , flf) 



+ A s c (a ->■ ag) 



a=q,q 

+(xi <-¥ x 2 )}dxidx 2 , 



G q / p {Xi,fif)Gq/ p (x 2 , (if) 



(34) 



where 



with 



Af{q -> qg) = Af{q g^) = C F (2 In <5 S + 3/2), 
^ c = Arin(4), 

G a (=q,q)/ P (x,fif) = 22 / —Gi3/p(x/y } n f )P a p(y) 



P a ?(y) = P aP {y) in(*c— 4) - 



(35) 
(36) 

(37) 



(38) 
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Finally, the NLO total cross section for pp — > ^Gkk in the MS factorization scheme is 

M 2 



a 



NLO 



dm S$-i 



M 2+5 m5 ' 1 ^ L ° 



(39) 



with 



NLO 



dx\dx 2 { 



G q /p(xi,fi f )Gg /p (x 2 ,fif) + (x 1 «-» x 2 ) 



coll 
m 



+ I (J ■>'] <!■>' J 

( a =g,0=q,q) 



G a / p (Xx, IIf)Gp/ p (x 2 , /If) + (Xx -H- X 2 ) 



&Z(a(3 -> -yG KK + 0) (40) 



Note that the above expression contains no singularities since 2A 2 + A 2 = and 1A\ + A\ + 
A f( q -> + Af (g -> = 0. 



IV. NUMERICAL RESULTS 



In this section we present numerical results for the total and the differential cross sections 
for 'jGkk associated production at the LHC In our numerical calculations the running QCD 
coupling constant a s (fi) is evaluated at three-loop order[40] and the CTEQ6.6M PDFs41| 
are used throughout. Only u and d flavor are activated since numerical calculations show 
that contributions from other flavor can be omitted. We take the LED parameters Mjj 
and 5 as input. Except for the scale uncertainty plot, both the renormalization and the 
factorization scale are fixed at pj,, which is the transverse momentum of the photon. Jet are 
defined by the following requirements: 



pf > 20GeV, 
\r] jet \ < 2.5. 



(41) 



Besides, the following cuts are assumed in our calculations [12j, |42 



_7 ^ m min 

Pt>Pt , 



\v\ < 2.4, 



miss — mm 

Pt Pt i 



&<P{l,PT iss ) > 2-5. 



(42) 



Here the default value of p™ m is 400GeV in the following calculations, as suggested in Ref. [12, 
421 ]. r\ is the pseudorapidity of the photon and p™ lss is the missing transverse momentum, 
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defined as 

p T , no jet in the final state, 



miss 

Pt 



Pt, with jet in the final state, 

where pj, is the transverse momentum of the graviton. We also require the photon to be iso- 
lated by requiring the separation of the photon and the radiated parton AR = a/ A</> 2 + Arj 2 
to be greater than 0.4. 

Moreover, it should be noted that the LED model is an effective low energy theory. 
Therefore we present two classes of numerical results to quantify the ultraviolet sensitivity: 
one with the truncation Tn^ G < M|>, vn 1 Q KK being the invariant mass of the graviton 
and the photon, while the other one is not truncated. As pointed out in Ref.[l2] if the two 
results significantly differ the contributions arising from regions above Md dominate, and 
the calculations are not under control but if they do not the LED model is viable. 

In Figj6] we show that it is reasonable to use the two cutoff phase space slicing method 
in our NLO QCD calculations; i.e., the dependence of the NLO QCD predictions on the 
arbitrary cutoffs S s and 5 C is indeed very weak, as was also found in Ref. 29[. While the 
Born cross sections and the virtual corrections are cutoff independent, both the soft and 
collinear contributions and the noncollinear contributions depend strongly on the cutoffs. 
However, the cutoff dependence in the two contributions (cr^ + a™ 11 and a^ c + <r^) nearly 
cancel each other , especially for the cutoff 5 S between 10 -4 and 10~ 3 , where the final results 
for <j NLO are almost entirely independent of the cutoffs. Therefore, we will take S s = 1CT 4 
in the numerical calculations below. Generally 5 C being 50 — 100 times smaller than 5 S 
is sufficient for accurate calculations to a few percent [29|, so we take 5 C = 5 s /50 in our 
calculations. 

Figure [7] shows the dependence of both the LO and the NLO total cross sections on 
the factorization scale(/i/) and the renormalization scale(/x r ) assuming Md = 2TeV, 5 = 4, 
and setting p™ n = 400GeV. When the scale \i varies from 0.2p T to 5p T , the LO total 
cross sections vary from 1.18 to 0.86fb, while the NLO total cross sections vary from 1.13 to 
0.91fb. Thus, the NLO corrections reduce the scale dependence, which makes the theoretical 
predictions somewhat more reliable. The conclusion is similar for 6 = 2, which is not shown 
here. 

In FigJH] and Figj9] we show the dependence of both the LO and the NLO total cross 
sections on M D , setting p mm = 400GeV, and assuming 5 = 2 and 4, respectively. As M D 
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increases the LO total cross sections decrease and the two results, with and without the 
truncation, approach each other. Also shown in Fig|8] and Figj9] are the K factors, defined 
as a /<j lo , which are around 1.3 ~ 1.5 for 5 = 2 and 1.1 ~ 1.3 for 5 = 4, respectively. 



We also give K factors for cases with jet veto[42j, where events with > lOOGeV are 
vetoed. In this case, the K factors are around 0.9 ~ 1 for 5 = 2 and 0.8 ~ 0.9 for 5 = 4, 
respectively, i.e., the NLO corrections reduce the LO results, which is due to the fact that 
the jet veto discards large positive contributions from real emission processes. 

In Fig{10] and Fig{TT] we show the dependence of both the LO and the NLO total cross 
sections on p™ n , for Md = 3TeV and 5 = 2 and 4, respectively. As p™ m increases the LO 
total cross sections decrease and the two results, with and without the truncation, differ 
increasingly. The K factors are about 1.3 ~ 1.4 for 5 = 2 and 1.2 for 5 = 4, respectively. 
When jet veto is considered, the K factors are around 0.95 for 5 = 2 and 0.84 ~ 0.9 for 
5 = 4, respectively. 

In Figs JT2lfT4l we display differential cross sections with truncation as functions of the 
missing transverse momentum, the transverse momentum and the pseudorapidity of photon, 
respectively. We find that the NLO QCD corrections always enhance the LO differential 
cross sections but do not significantly change the shapes of the LO differential cross sections. 

In conclusion, we have calculated the complete NLO QCD corrections to the inclusive 
total cross sections for jGkk associated production in the LED model at the LHC The 
NLO corrections generally enhance the total cross sections and the K factor is around 
1.3 ~ 1.5 for 5 = 2 and 1.1 ~ 1.3 for 5 = 4, respectively. When jet veto is considered, the 
NLO contributions reduce the LO results, the K factors are around 0.9 ~ 1 for 5 = 2 and 
0.8 ~ 0.9 for 5 = 4, respectively. We also compared the results with and without trunca- 
tion of tti^Gkk to quantify the ultraviolet sensitivity of the LED model. The NLO QCD 
corrections were found to reduce the dependence of the total cross sections on the renor- 
malization/factorization scale. We also calculated some important differential distributions 
for this process at the NLO, including the missing transverse momentum distribution, the 
transverse momentum distribution and the pseudorapidity distribution of photon. We found 
that the NLO corrections enhance these LO differential cross sections but do not appreciably 
change their shapes. 



15 



Acknowledgements 



This work is supported in part by the National Natural Science Foundation of China under 
grants No.10721063, No.10975004, and No.10635030, and the US Department of Energy, 
Division of High Energy Physics under Grant No. DE-FG02-91-ER4086. 



APPENDIX A 



In this appendix we give the related Feynman rules 



nidi 



q(k 2 )q(ki)G IMU : -i^llnih ~ k 2 ) v + lu(h - fc 2 ) M ], 

V a {k 1 )Vp{k 2 )G ^ -zf [h ■ k 2 C^ a p + D^ a p(ki, k 2 ) + E^^h, k 2 )}, 



In all the Feynman rules the particle momenta flow inward, gT a represents either 
g s T a if V is a gluon or eQf if V is a photon and 



D^apih, k 2 ) = r]^ u k 1/3 k 2a - [rj^k lu k 2a + ri^ a k lf) k 2u - r] a/3 k lfl k 2u + (fi -H- u)}, 
E^apih, k 2 ) = rj^vikiahp + k 2a k 2f3 + k la k 2f3 ) - [r] u/3 k lfl k la + T] ua k 2fl k 2 p + (ju -B- v)](43) 

APPENDIX B 

We collect the explicit expressions of the finite terms of the matrix element squared in 
this appendix. S in Eq. (fl2l) is given by 

k 2 G 2 1 

S =—^{ — [3C 2 t{m 2 - 4t)(m 4 - 2tm 2 + s 2 + t 2 )- 3C 5 (m 2 - 4t)(m 2 - t)(m 4 - 2tm 2 + s 2 
6 ut 

+ t 2 ) - 3Dlst(m 2 - 4t)(m 4 - 2tm 2 + s 2 + t 2 ) + 3Z^sw(2s 2 + 2ts + t 2 )(3m 2 - 4(s + t)) 
+ 3C 6 (2s 3 + 4ts 2 + 3t 2 s + t 3 )(3m 2 - 4(s + *)) + 3C 3 w(2s 2 + 2ts + t 2 )(4(s + 1) - 3m 2 ) 
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- 3C 4 (m 2 - s)(m 6 - Qtm 4 + (-5s 2 - Qts + 6t 2 )m 2 + 4s(2s 2 + 3ts + 3t 2 )) 

- 3C 1 ((s + 8t)m 6 - 6t(s + At)m 4 + (7s 3 + 18ts 2 + 30t 2 s + 32t 3 )m 2 - 4(2s 4 + 5ts 3 

+ 9tV + 8t 3 s + 4t 4 ))] + ^[2((3s + 14t)m 6 - 6(s 2 + 5ts + t 2 )m 4 + (3s 3 + 30ts 2 

- 16t 3 )m 2 + 2t(-7s 3 - 3ts 2 + 8t 2 s + 4t 3 )) + 18(-(s + 5t)m 6 + (s 2 + 8ts + lH 2 )m 4 

- (s 3 + 8ts 2 + 15t 2 s + 12t 3 )m 2 + t(5s 3 + llts 2 + 12t 2 s + 6t 3 ))] + 

stow-nL'-ins+ty ^ + 4ta ~ 2i2)m ' 2 " 2(3sS + 16ts2 + 17i2s 

+ 4t 3 )m 10 + (6s 4 + 22ts 3 + 56tV + 68t 3 s + 34t 4 )m 8 - 2(3s 5 - 5ts 4 - 41tV 

- 46tV + 5t 4 s + 18t 5 )m 6 + (3s 6 + 2ts 5 - llOtV - 300t 3 s 3 - 266tV - 60t 5 s 
+ 12t 6 )m 4 + 2st(-3s 5 + 14ts 4 + 86t 2 s 3 + 137t 3 s 2 + 84t 4 s + 16t 5 )m 2 

- 24s 2 t 3 (s + t) 3 ) log I — m 2 - 3s(m 2 - t) 2 (s + t) 2 (3m 10 - 6(s + 3t)m 8 + 6(s 2 + Ms 



+ 3t 2 )m 6 - 2s 2 (3s + llt)m 4 + s 2 (3s 2 + 42ts + 22t 2 )m 2 - 20s 3 t(s + *)) log 




- (m 2 - s)((m 2 - t)((s + t) 2 (21s + 34t)m 10 - (18s 4 + 196ts 3 + 483tV + 44H 3 s 

+ 136t 4 )m 8 + (s + t) 2 (15s 3 + 165ts 2 + 466t 2 s + 238t 3 )m 6 - (18s 6 + 124ts 5 + 587t 2 s 4 
+ 1473t 3 s 3 + 1786t 4 s 2 + 998t 5 s + 204t 6 )m 4 + t(s + t) 2 (49s 4 + 154ts 3 + 382tV 
+ 340t 3 s + 68t 4 )m 2 - 2st 2 (s + t) 3 (lls 2 + 34ts + 34t 2 ) - 3(m 2 - s)s(m 2 - t)t({2s 

+ 5t)m 6 - 6t(s + t)m A + 3(2s -t)(s + tfm 2 - 4(2s -t)(s + t) 3 ) log ) 

- 3st(s + t) 2 (3(4s + 3t)m 8 - (24s 2 + 57ts + 22t 2 )m 6 + (12s 3 + 75ts 2 + 74t 2 s 

+ 17t 3 )m 4 - t(27s 3 + 64ts 2 + 33t 2 s + 4t 3 )m 2 + 4st 2 (3s 2 + 4ts + t 2 )) log ( — V) )]}, 




(44) 
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where Q q is the electric charge of the initial (anti) quark, and 



7T 2 






3s' 






^ [ 




+ T» 


l^ 2 




\ 7T 2 



2(m 2 -s)' 

"■SRhShjT)^))-^ 

= 4h l0 8 2 f — ) + » 2 - 2L W " — ) " 2Li ^ - ^ )1. 
st \ s J s t 

Dl = -[-log 2 f— )+7r 2 -2Li 2 (l-— ) - 2Li 2 (l - — )]. (45) 
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FIG. 6: Dependence of the NLO total cross sections for the jGkk associated production at the 
LHC on the theoretical cutoff 5 S with 5 C = <5 s /50, assuming Mb = 3TeV, 5 = 4. Truncation 
mt n < is used here. 
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FIG. 7: Dependence of the NLO total cross sections for the ^Gkk associated production at the 
LHC on the factorization scale(/i/) and the renormalization scale(/U r ), assuming Mp = 2TeV, 
(5 = 4. Truncation mt n < M"k is used here. 
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FIG. 8: Dependence of the total cross section of the 'jGkk associated production at the LHC on 
Mn, assuming 5 = 2. 
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FIG. 9: Dependence of the total cross section of the 'jGkk associated production at the LHC on 
Mn, assuming 5 = 4. 
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FIG. 10: Dependence of the total cross section of the ^Gkk associated production at the LHC on 
p™ in , assuming M D = 3TeV, 5 = 2. 
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FIG. 11: Dependence of the total cross section of the ^Gkk associated production at the LHC on 
p™ in , assuming M D = 3TeV, 5 = 4. 
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FIG. 12: Dependence of the differential cross section of the ^Gkk associated production at the 
LHC on p™ ss , assuming Mb = 3TeV, 5 = 2 and 4, respectively. 
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FIG. 14: Dependence of the differential cross section of the ^Gkk associated production at the 
LHC on \r/\, assuming Md = 3TeV, 5 = 2 and 4, respectively. 
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